Expected Production of Strange Baryons and Antibaryons in Baryon-Poor QGP 
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In a dynamical model of QGP at RHIC we obtain the temporal evolution of strange phase space 
occupancy at conditions expected to occur in 100+lOOA GeV nuclear collisions. We show that the 
sudden QGP break up model developed to describe the SPS experimental results implies dominance 
of both baryon and antibaryon abundances by the strange baryon and antibaryon yields. 
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We explore the consequences of high strangeness abun- 
dance we expect to be present in the baryon-poor quark- 
gluon plasma (QGP) environment formed e.g., in the 
central rapidity region in Au-Au, maximum energy 
100+100^ GeV collisions at the relativistic heavy ion 
collider (RHIC) at the Brookhaven National Laboratory, 
Upton, New York. About 10-20% of hadrons produced 
in these reactions will be strange, and since mesons dom- 
inate hadron abundance, there is much more strangeness 
than baryon number. During the break-up of the color 
charge deconfined QGP phase there is considerable ad- 
vantage for strangeness to stick to baryons given that 
the energy balance for the same flavor content favors 
production of strange baryons over kaons, {E{K + -k) < 
i?(N+K)) . When QGP is formed, we therefore expect 
to find hyperon dominance of baryon distribution i.e. 
most baryons and antibaryons produced at RHIC will 
be strange. A similar argument could be made for re- 
actions leading to the confined phase, the main differ- 
ence arises from the observation that the required high 
abundance of strangeness per participating nucleon can 
be produced in the deconfined QGP This qualita- 

tive argument will be quantitatively elaborated here, in 
view of the considerable effort that has been committed 
by the STAR collaboration at RHIC to enhance the capa- 
bility to measure multi-strange (antijbaryon production 
using a silicon strip detector (SSD) g. 

In the first part of this report we show that the chem- 
ical strangeness flavor abundance equilibrium is estab- 
lished at the time of QGP breakup by the dominant pro- 
cess which is gluon fusion, GG — > ss. In QGP this re- 
action overcomes the current quark mass threshold 2tos 
for strangeness formation for temperature T ~ mg > 
200 MeV. In the second part of this report we use the 
computed chemical condition of strangeness, and employ 
the knowledge gained in our analysis of the SPS results 
1^, to obtain the strange baryon and antibaryon abun- 
dances expected at RHIC. 

In some key aspects the methods we employ differ from 
those obtained in other studies of chemical equilibration 
of quark flavor for RHIC conditions j7|H] . We study the 
dynamics of the phase space occupancy rather than par- 



ticle density, and we eliminate most of the dynamical 
flow effects by considering entropy conserving evolution. 
Moreover, we use running QCD parameters (both cou- 
pling and strange quark mass) to describe strangeness 
production, with strong coupling constant as as deter- 
mined at the Mzo energy scale. We will make two as- 
sumptions of relevance for the results we obtain: 

a) the kinetic (momentum distribution) equilibrium is 
reached faster than the chemical (abundance) equilib- 
rium 0,|ll|; 

b) gluons equilibrate chemically significantly faster than 
strangeness | |T^ . 

The first assumption allows to study only the chemical 
abundances, rather than the full momentum distribution, 
which simplifies greatly the structure of the master equa- 
tions; the second assumption allows to focus after an ini- 
tial time To has passed on the evolution of strangeness 
population: tq is the time required for the development 
to near chemical equilibrium of the gluon population. 

We now formulate the dynamical equation for the evo- 
lution of the phase space occupancy 7s of strange quarks 
in the expanding QGP: the phase space distribution /s 
can be characterized by a local temperature T{x, t) of a 
(Boltzmann) equilibrium distribution , with normal- 
ization set by a phase space occupancy factor: 



/,(p,f;i))c.7,(T)/-(p;r). 



(1) 



Eq.(0) invokes in the momentum independence of 7s 
our first assumption. More generally, the factor 7i, i = 
g, q, s, c allows a local density of gluons, light quarks, 
strange quarks and charmed quarks, respectively not to 
be determined by the local momentum shape, but to 
evolve independently. 

With variables {t, x) referring to an observer in the lab- 
oratory frame, the chemical evolution can be described by 
the strange quark current non-conservation arising from 
strange quark pair production described by a Boltzmann 
collision term: 
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The factor 1/2 avoids double counting of gluon pairs. 
The impHcit sums over spin, color and any other discreet 
quantum numbers are combined in the particle density 
p = ^ J <Pp /, and we have also introduced the mo- 
mentum averaged production/annihilation thermal reac- 
tivities: 



Jd^PlJd^P2f{pi,T)f{p2,T) 



(3) 



f(j)i,T) are the relativistic Boltzmann/Jiittncr distribu- 
tions of two colliding particles z = 1, 2 of momentum pi. 

The current conservation can also be written with ref- 
erence to the individual particle dynamics ]l3t : con- 
sider ps as the inverse of the small volume available to 
each particle. Such a volume is defined in the local 
frame of reference for which the local flow vector van- 
ishes 'ij{x,t)\^Q^^^ = 0. The considered volume 6Vi being 
occupied by small number of particles 6N {e.g., SN = 1), 
we have: 



(4) 



The left hand side (LHS) of Eq. (g) can be now written 
as: 

dps dvps _ 1 dSNs _ dps 1 dSVi 

dt dx ~ 5Vi dt ~ dt "^"SVi dt ' ^ ' 

Since SN and SVidt are L(orentz)-invariant, the actual 
choice of the frame of reference in which the right hand 
side (RHS) of Eq. (|) is studied is irrelevant and we drop 
henceforth the subscript I. 

We can further adapt Eq. to the dynamics we pur- 
sue: we introduce p^{T) as the (local) chemical equilib- 
rium abundance of strange quarks, thus p — 7sP^. We 
evaluate the equilibrium abundance SN^ = 6Vp^{T) 
integrating the Boltzmann distribution: 



SN,"" = [5VT^]^z^K2{z), 



T 



(6) 



where Ki, is the modified Bessel function of order we 
will below use: d[z'^ K,y{z)]/dz = —z^K^-i . The first 
factor on the RHS in Eq. (g) is a constant in time should 
the evolution of matter after the initial pre-thermal time 
period Tq be entropy conserving |p^ , and thus 5VT^ = 
JVqTq =Const. . We now substitute in Eq. (|^) and obtain 

dp^ , dvps f d-y^ ^ -fs_ Ki{z) \ 

dt^ dx " \dT ^ T K2{z)) ' ^ ' 

where T ~ dT/dt. Note that in Eq. only a part of the 
usual flow-dilution term is left, since we implemented the 
adiabatic volume expansion, and study the evolution of 
the phase space occupancy in lieu of particle density. The 
dynamics of the local temperature is the only quantity 
we need to model. 



We now return to study the collision terms seen on the 
RHS of Eq. (||). A related quantity is the (L- invariant) 
production rate ^12^34 particles per unit time and 
space, defined usually with respect to chemically equili- 
brated distributions: 
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The factor 1/(1 + 61,2) is introduced to compensate 
double-counting of identical particle pairs. In terms of 
the L-invariant A , Eq. (||) takes the form: 

\dT T K2{z) J '9^ ' 
+7girhg{T)A'i'^''- 7,(T)7g(r)(A-"-99+ ^"^-??). (9) 

Only weak interactions convert quark flavors, thus, on 
hadronic time scale, we have 7s,ij(t) = 7s,q(T). Moreover, 
detailed balance, arising from the time reversal symme- 
try of the microscopic reactions, assures that the invari- 
ant rates for forward/backward reactions are the same, 
specifically 
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and thus: 
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When all 7^ 1, the Boltzmann collision term vanishes, 
we have reached equilibrium. 

As discussed, the gluon chemical equilibrium is 
thought to be reached at high temperatures well before 
the strangeness equilibrates chemically, and thus we as- 
sume this in what follows, and the initial conditions we 
will study refer to the time at which gluons are chemically 
equilibrated. Setting Ag = 1 (and without a significant 
further consequence for what follows, since gluons dom- 
inate the production rate, also \q — 1) we obtain after 
a straightforward manipulation the dynamical equation 
describing the evolution of the local phase space occu- 
pancy of strangeness: 
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Here, we defined the relaxation time Tg of chemical 
(strangeness) equilibration as the ratio of the equilibrium 
density that is being approached, with the rate at which 
this occurs: 
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The factor 1/2 is introduced by convention in order for 
the quantity Tg to describe the exponential approach to 
equiUbrium. 

Eq. ( [T^ ) is our final analytical result describing the evo- 
lution of phase space occupancy. Since one generally ex- 
pects that 7s ^ 1 in a monotonic fashion as function of 
time, it is important to appreciate that this equation al- 
lows 7s > 1: when T drops below TOs, and 1/ts becomes 
small, the dilution term (2nd term on LHS) in Eq. ( p^ 
dominates the evolution of 7s . In simple terms, the high 
abundance of strangeness produced at high temperature 
over-populates the available phase space at lower tem- 
perature, when the equilibration rate cannot keep up 
with the expansion cooling. This behavior of 7s has been 
shown in Fig. 2 of Ref jl^ for the SPS conditions with fast 
transverse expansion. Since we assume that the dynam- 
ics of transverse expansion of QGP is similar at RHIC 
as at SPS, we will obtain a rather similar behavior for 
7s. We note that yet a faster transverse expansion than 
considered here could enhance the chemical strangeness 
anomaly. 

Ts (T) , Eq. ( |T3| ) , has been evaluated using pQCD cross 
section and employing NLO (next to leading order) run- 
ning of both the strange quark mass and QCD-coupling 
constant as jl^. We believe that this method produces 
a result for as that can be trusted down to 1 GeV en- 
ergy scale which is here relevant. We employ results ob- 
tained with as{Mzo) = 0.118 and ms(lGeV) = 220 MeV, 
a somewhat conservative (high) choice for TOs, which 
should under-predict strangeness production. There is 
some systematic uncertainty due to the appearance of the 
strange quark mass as a fixed rather than running value 
in both, the chemical equilibrium density in Eq. (p^), 
and in the dilution term in Eq. (p^). We use the value 
ms(lGeV), with the energy scale chosen to correspond 
to typical interaction scale in the QGP. 

Numerical study of Eq. (|l^ ) becomes possible as soon 
as we define the temporal evolution of the temperature 
for RHIC conditions. We expect that a global cylindri- 
cal expansion should describe the dynamics: aside of the 
longitudinal flow, we allow the cylinder surface to expand 
given the internal thermal pressure. SPS experience sug- 
gests that the transverse matter fiow will not exceed the 
sound velocity of relativistic matter v±^ ~ c/%/3- We re- 
call that for pure longitudinal expansion local entropy 
density scales as S" cx T'^ c>c 1/t, [^. It is likely that 
the transverse flow of matter will accelerate the drop in 
entropy density. We thus consider the following temporal 
evolution function of the temperature: 



T{t) = To 



-,1/3 



[l+T 2c/d){l+T v^/R^)^ 



(14) 



We take the thickness of the initial collision region at 
To = 0.5 GeV to be d{To = 0.5)/2 = 0.75 fm, and the 
transverse dimension in nearly central Au-Au collisions 



to be R± = 4.5 fm. The time at which thermal initial con- 
ditions are reached is assumed to be tq = Ifm/c. When 
we vary Tq, the temperature at which the gluon equilib- 
rium is reached, we also scale the longitudinal dimension 
according to: 



d(To) = (0.5GeV/ro)^1.5fm. 



(15) 



This assures that when comparing the different evolu- 
tions of 7s we are looking at an initial system that has 
the same entropy content. The reason we vary the ini- 
tial temperature Tq down to 300 MeV, maintaining the 
initial entropy content is to understand how the assump- 
tion about the chemical equilibrium of gluons, reached 
by definition at Tq, impacts our result. 

The numerical integration of Eq. (|l^) is started at tq, 
and a range of initial temperatures 300 < Tq < 600, 
varying in steps of 50 MeV. The high limit of the tem- 
perature we explore exceeds somewhat the "hot glue sce- 
nario" ||l^, while the lower limit of Tq corresponds to 
the more conservative estimates of possible initial con- 
ditions |l^ . Since the initial p-p collisions also produce 
strangeness, we take as an estimate of initial abundance 
a common initial value 7s (To) = 0.2. The time evolution 
in the plasma phase is followed up to the break-up of 
QGP. This condition we establish in view of our analysis 
of the SPS results. We recall that SPS-analysis showed 
that the system dependent baryon and antibaryon m±- 
slopes of particle spectra are result of differences in collec- 
tive flow in the deconflned QGP source at freeze-out 0. 
There is a universality of physical properties of hadron 
chemical freeze-out between different SPS systems, and in 
our analysis a practical coincidence of the kinetic freeze- 
out conditions with the chemical freeze-out. We thus 
expect extrapolating the phase boundary curve to the 
small baryochemical potentials that the QGP break-up 
temperature T^^^ ~ 145 ± 5 MeV will see just a minor 
upward change to the value T^^hic ^ _^ 5 MeV. 

With the freeze-out condition fixed, one would think 
that the major uncertainty in our approach comes from 
the initial gluon equilibration temperature Tq, and we 
now study how different values of Tq inffuence the fi- 
nal state phase space occupancy. We integrate numeri- 
cally Eq. ( [12I ) and present 7s as function of both time t in 
Fig.j^a, and temperature T in Fig. ^d, up to the expected 
QGP breakup at T^hic 159 ± 5 MeV. We see that: 

• widely different initial conditions (with similar initial 
entropy content) lead to rather similar chemical condi- 
tions at chemical freeze-out of strangeness, 

• despite a series of conservative assumptions we find not 
only that strangeness equilibrates, but indeed that the 
dilution effect allows an overpopulation of the strange 
quark phase space. For a wide range of initial conditions 
we obtain a narrow band 1.18 > 7s(7/) > 0.95 . We will 
in the following study of strange baryon and antibaryon 
abundances adopt what we believe to be the most likely 
value 7s(T/) = 1.15. 
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FIG. 1. Evolution of QGP-phase strangeness phase space 
occupancy 7s a) as function of time and b) as function of 
temperature, see text for details. 

We now consider how this relatively large value of 7s, 
characteristic for the underlying QGP formation and evo- 
lution of strangeness, impacts the strange baryon and 
anti-baryon observable emerging in hadronization. Re- 
membering that major changes compared to SPS should 
occur in rapidity spectra of mesons, baryons and an- 
tibaryons, we will apply the same hadronization model 
that worked in the analysis of the SPS data |^ : 

1) the QGP freeze-out/break-up occurs without a signif- 
icant (transient) hadronic gas epoch; 

2) the deconfined QGP state evaporates over a few fm/c, 
during which time it remains near to the freeze-out tem- 
perature, with energy lost due to particle evaporation 
and work done against the vacuum balanced by the in- 
ternal energy flows. This reaction picture can be falsified 
easily, since we expect, based and compared to the Pb- 
Pb 158A GeV results: 



a) shape identity of RHIC ■m± and y spectra of an- 
tibaryons p , A , S , since in our approach there is no 
difference in their production mechanism, and the form 
of the spectra is determined in a similar way by the local 
temperature and flow velocity vector; 

b) the m_L-slopes of these antibaryons should be very sim- 
ilar to the result we have from Pb~Pb 158 A GeV since 
only a slight increase in the freeze-out temperature oc- 
curs, and no increase in collective transverse flow is ex- 
pected. 

The abundances of particles produced from QGP 
within this sudden freeze-out model are controlled by 
several further chemical parameters: the light quark fu- 
gacity 1 < Xq < 1.1, value is hmited by the expected 
small ratio between baryons and mesons (baryon-poor 
plasma) when the energy per baryon is above 100 GeV, 
strangeness fugacity As — 1 which value for locally neu- 
tral plasma assures that (s — s) = 0; the light quark 
phase space occupancy 7^ ~ 1.5, overabundance value 
due to gluon fragmentation. Given these narrow ranges 
of chemical parameters and the freeze-out temperature 
Tf = 150 MeV, we compute the expected particle pro- 
duction at break-up. In general we cannot expect that 
the absolute numbers of particles we find are correct, as 
we have not modeled the important effect of flow in the 
laboratory frame of reference. However, ratios of hadrons 
subject to similar flow effects (compatible hadrons) can 
be independent of the detailed final state dynamics, as 
the results seen at SPS suggest Q, and we will look at 
such ratios more closely. 

Taking 7, = 1.25, 1.5, 1.6 we choose the value of Xq, 
see the header of table ^, for which the energy per baryon 
{E/b) is similar to the collision condition (100 GeV), 
which leads here to the range Xq = 1.03 ± 0.005. We 
evaluate for these examples aside o{ E/b, the strangeness 
per baryon s/b and entropy per baryon S/b as shown 
in the top section of the table |. We do not enforce 
(s — s) = exactly, but since baryon asymmetry is small, 
strangeness is balanced to better than 2% in the param- 
eter range considered. In the bottom portion of table 
I we present the compatible particle abundance ratios, 
computed according to the procedure developed in 
We have presented aside of the baryon and antibaryon 
relative yields also the relative kaon yield, which is also 
well determined within our approach. 

The meaning of these results can be better appreciated 
when we assume in an example the central rapidity den- 
sity of protons is dp/ dylcentrai = 25. In table || we present 
the resulting (anti)baryon abundances. We see that the 
net baryon density db/dy ~ 15 ± 2, there is baryon num- 
ber transparency. We see that (anti)hyperons are indeed 
more abundant than non-strange (anti)baryons. It is im- 
portant when quoting results from table ^ to recall that: 

1) we have chosen arbitrarily the overall normalization in 
table only particle ratios were computed, and 

2) the rapidity baryon density relation to rapidity pro- 



4 



ton density is a consequence of the assumed value of A,, 
which we chose to get E/b ~ 100 GeV per participant. 

However, we firmly believe that our key result seen in 
table H , the hyperon-dominance of the baryon yields, 
does not depend on these hypothesis. Indeed, we have 
explored another set of parameters in our first and pre- 
liminary report on this matter , finding the same pri- 
mary conclusion explained in the introduction. Another 
related notable result, seen in table |[ is that strangeness 
yield per participant is 13-23 times greater than seen 
at present at SPS energies, where we have 0.75 strange 
quark pairs per baryon. As seen in table || the baryon 
rapidity density is in our examples similar to the proton 
rapidity density. 

In summary, we have shown that one can expect 
strangeness chemical equilibration in nuclear collisions 
at RHIC if the deconfined QGP is formed, with a 
probable overpopulation effect associated with the early 
strangeness abundance freeze-out before hadronization. 
We studied the physical conditions at QGP breakup and 
have shown that (anti)hyperons dominate (anti)baryon 
abundance. 
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TABLE I. For 7s = 1.15, As = 1 and 7,, Xq as shown: 
Top portion: strangeness per baryon s/b, energy per baryon 
i?/6[GeV] and entropy per baryon S/b. Bottom portion: sam- 
ple of hadron ratios expected at RHIC. 
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